H-function extension of the NBD: further applications 
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The H-function extension of the Negative Binomial Distribution is investigated for scaling expo- 
nents fi < 0. Its analytic form is derived via a convolution property of the H-function. Applications 
are provided using multihadron and galaxy count data for P n . 
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In a series of recent papers the present author developed a generalization of the Negative Binomial Distribution [1-3] . 
P n was obtained by means of the Poisson transform [4] 
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of the continuous probability density f(x) which asymptotically yields the KNO scaling form of P n (provided that 
f(x) belongs to the family of scale parameter generating density functions). Choosing f(x) to be the generalized 
gamma distribution [5] 
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one obtains the above cited generalization of the NBD. In Eqs. (2-3) we have parameters k > (shape parameter), 
A > (scale parameter), and /1 ^ (scaling exponent). In the limit fi — > and k — > 00 the generalized gamma density 
converges to the log-normal law. The k = 1, fi > special case of Eq. (2) provides the Weibull distribution whereas 
the /x = 1 case is, of course, the ordinary gamma distribution whose Poisson transform yields the NBD. 

In refs. [1-3] the Poisson transform of the generalized gamma density Eq. (2) was investigated only for scaling 
exponents /1 > 0. The analytic form of P n was shown to be expressible in terms of Fox's H-function (see refs. [1,3] 
for a summary and ref. [6] for details). Therefore we will call this particular extension of the Negative Binomial 
Distribution as HNBD for short. It takes the form 
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where H(-) denotes the H-function of Fox and 9 — A 1 /' 1 can be expressed in terms of the average multiplicity (n) 
according to 9 — (n)T(k)/T(k + l//x). The pure NBD given by Eq. (4) is the fj, = 1 marginal case of Eq. (5) for 
(ri) < k and of Eq. (6) for (n) > k [3]. The /1 — * 0, k — > 00 Poisson transformed log-normal limit of the HNBD lacks 
a representation in terms of known functions. 

The main goal of the present Letter is to complete the derivation of the Poisson transformed generalized gamma 
distribution by extending its validity to negative values of the scaling exponent. For /1 < we shall consider the 
Poisson transform of the probability density 
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f(x) = W) xkx ~' tk ~ 1 exp (~ Aar ^ 

with /i > 0. In the derivation of P„ two particular cases of H(a;) will be utilized, 
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the latter provides xf(x) with = A 1 /^ as before. We shall also make use of the following two properties of Y\(x): 
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with parameters 
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For the conditions of validity of various properties of H(x), see ref. [6]. Absorbing x n into the lhs. of Eq. (9) via 
identity (10) one obtains for P n defined by Eq. (1) with f(x) given by Eq. (7) the following integral: 
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Comparison with the convolution property Eq. (11) yields 
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for the Poisson transform of the probability density Eq. (7). Above, is expressible in terms of (n) according to 
9 = (n)T(k)/T(k — 1/fi) if k > This can be seen from Eq. (3) recalling that (x q ) = {n(n — 1) . . . (n — q + 1)) for 
the Poisson transform Eq. (1). Accordingly, the factorial moments of Eq. (15) may diverge. The generating function 
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can also be derived via the convolution property Eq. (11), it turns out to be 
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Changing u to —it in Eq. (17) the characteristic function (p(t) = J °° e ltx f(x) dx is obtained for Eq. (7). The generating 
function of the HNBD corresponding to Eqs. (5-6) was determined in [1]. 
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With Eqs. (15,17) we have completed the derivation of the Poisson transformed generalized gamma distribution. 
According to a theorem of Bondesson, for > 1 the characteristic function of Eq. (2) is an entire analytic function of 
finite order and therefore it must have complex zeroes which is not permitted for infinitely divisible entire characteristic 
functions [7]. Since the infinite divisibility of f(x) is preserved by P n for Eq. (1) one can deduce that this important 
feature holds for the HNBD if < < 1. In rcf. [1] it has already been demonstrated that the factorial cumulant 
moments of the HNBD exhibit nontrivial (not alternating) sign-changing oscillations for fx > 1 in accordance with 
the violation of infinite divisibility of P n . The effect is illustrated in Fig. 1 for shape parameter k = 1 (Wcibull case) 
which describes the inelastic pp and deep-inelatic e + p multipilicity data very well [2,3]. Bondesson's theorem tells us 
that sign-changing oscillations of the factorial cumulants may arise also for ji < — 1. But some preliminary results 
indicate that this is not the case and the factorial cumulants of the HNBD, when exist, remain positive for /j, < — 1. 

Analysing the experimental data for P n in different collision processes it was found that the HNBD shape with 
H < is, although rare, not completely absent. The best example is the full phase-space TASSO data in e + e~ 
annihilations at yfs = 34.8 GeV [8]. Most of the theoretical models struggle in the description of this high statistics 
data set. One of the few exceptions is ref. [9] reporting x 2 /d.o.f. = 9.1/14. It corresponds to the NBD fitted in two 
matched domains of multiplicity n with 3 fit parameters altogether. The HNBD analysis was carried out with fixed 
shape parameter k = 1 (details of the fitting procedure can be found in [1,3]). This particular case of the HNBD 
produces x 2 / d - oi - = 7.9/15 with parameters (n) = 13.595 ± 0.050 and /x = -10.053 ± 0.388. The fit is illustrated in 
Fig. 2. At Vs = 22 GeV we obtained x 2 /d.o.f. = 4.9/11, (n) = 11.313 ± 0.096 and fi = -12.745 ± 1.449. Despite 
of the relatively small statistics available at lowest PETRA energy y/s = 14 GeV, the HNBD fit to the TASSO data 
fails with x 2 /d.o.f. = 20.4/10. At top PETRA energy yfs = 43.6 GeV the fj, — » 0, k — > oo log-normal limit of the 
HNBD produces the best quality fit (x 2 /d.o.f. = 6.4/16) similarly to previous results obtained at the Z° peak [1,3]. 
This change in the parametrization of the best-fit HNBD indicates a slight narrowing of the KNO functions with 
increasing s in the energy range investigated. It is worth mentioning that the value of \x was found in the fitting 
procedures to be extremely sensitive to tiny details of the experimental P n . Since the e + e~ multiplicity data exhibit 
approximate KNO scaling and dominantly follow the \i — Poisson transformed log-normal shape, the \i < type 
behaviour with divergent higher-order factorial moments should not be taken too seriously. 

Let us now turn our attention to a completely different application of the HNBD. It is widely known that the 
/i = 1 special case, the pure NBD, yields satisfactory results for galaxy count distributions. In ref. [10] the Zwicky 
catalogue of galaxy clusters was investigated. The probability of finding n galaxies in a Zwicky cluster ((n) ~ 100) 
was found to follow Eq. (2) on a KNO plot with \i = 1 and k <~ 5. We have performed a similar analysis using the 
Lick galaxy catalogue which is one of the most successful attempts to map the extragalactic sky. It is known also as 
Shane- Wirtanen survey after the name of the authors of the original catalogue [11]. The galaxy counts are contained 
in 1246 photographic plates covering <~ 70 % of the sky, essentially the same region as the Zwicky sample. Each plate 
covers 6° x 6° and consists of a 36 x 36 array of counts in 10' x 10' cells. The catalogue contains ~ 1.25 x 10 6 galaxies 
out of which cca. 34 % are double counts because each plate overlaps its neighbour by at least 1°. Most of the double 
counts can be eliminated by using only the center 5° x 5° of each plate. In ref. [12] several multiplicative correction 
factors were determined to the raw counts of Shane and Wirtanen. The gray-scale map of the corrected galaxy counts 
presented in [12] shows an impressive network of galaxies with lots of clumps and filaments. An early computer model 
universe designed to reproduce this fine structure can be found in [13] where the frequency distribution of galaxy 
counts in the 10' x 10' cells is also displayed. 

The count distribution of galaxies analysed here was extracted from the Lick catalogue in ref. [14]. The plate 
selection criteria, applied in statistical analyses of the catalogue to exclude plates near the galactic plane, resulted 
467 plates with 420300 cells in total. The average count in cells is (n) <~ 1.4. Despite of the small mean of counts, 
cells with n > 20 galaxies are not rare in the Lick survey and the count distribution has a very long tail extending 
well over n/(n) = 10. This is a sure sign of highly non-Poissonian behaviour of fluctuations. In [13] a comparison 
of the Shane- Wirtanen map and a completely random distribution of galaxies is shown — the difference between the 
visual appearance of the two patterns is huge. Another way of demonstrating this strong departure is fitting the Lick 
galaxy count distribution in the 10' x 10' cells by the fi = 1, k — > oo Poisson limit of the HNBD. The fit is displayed 
in the top left plot in Fig. 3. The best-fit Negative Binomial with k ~ 2 improves on the Poissonian but the deviation 
remains substantial, see the top right plot of the same figure. We mention that the analysis of the Zwicky catalogue 
in [10] also signals a similar problem with the u = 1 NBD/gamma model since the good agreement mentioned above 
was obtained after truncating the high-n tail of the data. In ref. [3] it was shown that long tailed distributions such 
as the P n in pp collisions at \fs = 900 GeV can be successfully described by the u — > 0, k — ► oo limit of the HNBD. 
Fitting the Poisson transformed log- normal law to our galaxy data one obtains further improvement, see the bottom 
left plot in Fig. 3, but the agreement in the high-n tail is, again, unsatisfactory. The main difficulty of reproducing the 
Lick counts lies in the fact that the distribution decays according to an inverse power-law. The bottom right plot of 
Fig. 3 illustrates a HNBD fit with negative scaling exponent fj,. As is seen this attempt finally succeeds in describing 
the high-n tail too. The best-fit parameters are u = -0.571 ± 0.070 and k = 8.512 ± 0.447. 
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Summarizing our results, we have investigated the Poisson transform of the generalized gamma distribution Eq. (2) 
for negative values of the parameter \i. The analytic form of P n was derived by considering the Poisson transform 
of the probability density Eq. (7) with fj, > 0. Making use of the convolution property Eq. (11) of the H-function 
we have obtained Eqs. (15) and (17) for the analytic form of P n and Q(u). Our result completes the specification 
of the Poisson transformed generalized gamma distribution introduced in [1] and developed in [2,3]. Its application 
to multihadron and galaxy count data revealed /i < type behaviour in both cases. Nevertheless, it seems to be 
unlikely that the HNBD shape with negative /i will be frequently encountered for the full phase-space multiplicity 
distributions in particle and nuclear collisions. This feature is more probable for galaxy count data which often decay 
according to an inverse power-law reflecting the scale-invariant distribution of galaxies in the universe. Concerning 
multiparticle dynamics, similar behaviour of P n should be searched in restricted domains of phase-space. 
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FIG. 1. Sign-changing oscillations of the factorial cumulant-to-moment ratios H q [15] for the 
k = 1 special case of HNBD. Considering slices with fixed rank q, the neighbouring bumps along 
the /i-axis correspond to H q having opposite sign. For /i < 1 the moment ratios are always positive. 
The /x-scale is logarithmic and for clarity log \H q \ is displayed only for odd ranks q. 




FIG. 2. The best-fit HNBD to the full phase-space TASSO data at = 34.8 GeV [8]. The fit 
parameters are: k = 1 fixed, (n) = 13.595 ± 0.050 and fi = -10.053 ± 0.388. 
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FIG. 3. Various special and limiting cases of HNBD fitted to the Lick galaxy count distribution. 
Top left: best-fit Poisson (// = 1, k — » oo); top right: best-fit Negative Binomial (/x = 1, k ~ 2); 
bottom left: best-fit Poisson transformed log-normal (// — » 0, k — » oo); bottom right: best-fit HNBD 
with parameters fx = -0.571 ±0.070 and k = 8.512 ±0.447. The average count in cells is (n) ~ 1.4. 
The errorbars represent Poisson errors. 
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